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Abstract. The local magnetization in the one-dimensional random-field Ising 
model is essentially the sum of two effective fields with multifractal probability 
measure. The probability measure of the local magnetization is thus the 
convolution of two multifractals. In this paper we prove relations between the 
multifractal properties of two measures and the multifractal properties of their 
convolution. The pointwise dimension at the boundary of the support of the 
convolution is the sum of the pointwise dimensions at the boundary of the support 
of the convoluted measures and the generalized box dimensions of the convolution 
are bounded from above by the sum of the generalized box dimensions of the 
convoluted measures. The generalized box dimensions of the convolution of 
Cantor sets with weights can be calculated analytically for certain parameter 
ranges and illustrate effects we also encounter in the case of the measure of the 
local magnetization. Returning to the study of this measure we apply the general 
inequalities and present numerical approximations of the D 9 -spectrum. For the 
first time we are able to obtain results on multifractal properties of a physical 
quantity in the one-dimensional random-field Ising model which in principle could 
be measured experimentally. The numerically generated probability densities 
for the local magnetization show impressively the gradual transition from a 
monomodal to a bimodal distribution for growing random field strength h. 



PACS numbers: 05.45.Df, 05.50.+q, 75.10.Nr, 05.70.Fh 

1. Introduction 

Multifractal measures appear in a variety of contexts. The one-dimensional random- 
field Ising models [l]-[20] and random exchange [21]-[23] Ising models as well as other 
one-dimensional disordered systems [24]- [27], Bernoulli convolutions [28] and even 
learning in neural networks [29]- [31] are prominent examples. The use of a reduction 
scheme for the partition function of the one-dimensional random-field Ising model first 
introduced by Rujan [10] results in the partition function of a one spin system in an 
effective field (one-sided reduction, spin at the boundary) [4]- [16] or in two effective 
fields (two-sided reduction, spin in the bulk) [17] which is the appropriate point of 
view when investigating the local magnetization. The effective fields are governed by 
an iteration thus giving rise to a random iterated function system which is known to 
have a unique invariant measure [32]. The invariant measure of the effective fields is 
typically a multifractal [1, 4, 6]. 

In a recent publication [9] phase transitions in the Z3 g -spectrum of this invariant 
measure of the effective field were investigated (cf also [11, 12]) and tight bounds on the 
D q based on the pointwise dimension at specific points generalizing results in [13] were 
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formulated. The combination of both allows a more or less complete understanding of 
the -D g -spectrum of the invariant measure of the effective field by exclusively analytical 
methods. Naturally the question arises whether these results are relevant for physical 
quantities such as correlation functions or the local magnetization which in principle 
are experimentally accessible. 

The local magnetization can be expressed as a function of the effective fields 
[4]- [6] and we show in this paper that a considerable amount of the knowledge of the 
multifractal properties of the invariant measure of the effective field can be transferred 
to the measure of the local magnetization. Being essentially the sum of two effective 
fields with multifractal probability measure the local magnetization has a probability 
measure which is the convolution of two multifractal measures. We therefore first 
prove general relations between the multifractal properties of two measures and the 
multifractal properties of their convolution which then can be applied to the random- 
field Ising chain. As the convolution of measures is the composition of constructing 
the product measure and projecting it in a certain way, the work on projections of 
multifractal (product) measures [33, 34] is related to our subject. Whereas these 
papers focus mainly on properties of projected measures with respect to typical 
projections we are here concerned with a given projection leading to the convolution. 
This special case need not necessarily have the properties of a generic projection. 
There is also related work on the superposition of multifractals [35]- [37] and some 
remarks on the convolution of multifractals in [38]. 

In addition to the mathematical results we also calculate the measure of the 
local magnetization and its multifractal spectrum numerically, random-field Ising 
systems can be realized as dilute antiferromagncts in uniform magnetic fields [39] 
and the local magnetization can in principle be measured by neutron scattering or 
MoBbauer spectroscopy. The probability distribution of the local magnetization with 
respect to the disorder therefore should be experimentally accessible and could be 
compared to our numerical results presented in figure 9. Especially the gradual 
transition from a strongly peaked monomodal distribution to a strongly peaked 
bimodal distribution observed numerically should clearly be visible. Depending on 
the quality of the measurement it is even feasible to calculate the Z? g -spectrum of the 
obtained probability distribution by the box methods described in section 5 and to 
compare to the results presented here. This should at least reproduce the general form 
of the multifractal spectrum shown in figure 7. 

The paper is organized as follows. After recalling the model and the reduction 
scheme in section 2 we prove general bounds on the D g -spectrum of the convolution 
of two measures and relations between certain pointwise dimensions in section 3. The 
results are applied to the situation of the local magnetization in the random-field Ising 
chain. We then explicitly calculate some D g -spectra in the simplified situation of the 
convolution of equal-scale Cantor sets with weights in section 4. In section 5 we present 
numerical results for the D^-spectrum of the measure of the local magnetization and 
in the concluding section 6 we summarize our results and draw some conclusions. 
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2. The model 

In the following we consider the one-dimensional random-field Ising model with 
quenched disorder which for N spins has the Hamiltonian 

6-1 b 

H N {{s} N ) = -J^2siS i+1 -^2hiSi, (2.1) 

i— a i—a 

with a < < b and b — a + 1 = N. s, denotes the classical spins at site i of the chain 
taking values ±1, J is the coupling strength between spins and hi is the random field 
at site i. The random fields are independent identically distributed random variables 
with probability density 

p{hi) = ^S(hi-h) + ^5(hi + h), heR+. (2.2) 

In former work we used a reformulation of the canonical partition function Zn = 
J2{ s } n ex P(~ (3Hn({s}n)) to the partition function of the spin s a at the left-hand 

boundary of the chain in an effective field which was first introduced by Rujan 
[10], 



N 

""=±1 



£ + 5]B(ij'°)]), (2.3) 



f' = <) + ^ 45 = 0, (2.4) 



with 



Pn{x) 



A(x) = (2/3)- 1 ln(cosh/3(x + J)/ cosh j3(x- J)), (2.5) 
B{x) = (2/3) -1 ln(4cosh/3(x + J) cosh/3(x - J)). (2.6) 

When viewing (2.4) as a random iterated function system (RIFS) we will also write 
x n instead of for the effective field after n = N — i + 1 iterations of (2.4). 

The effective fields random variables on the random-field probability space 

and we write p n (x) for their induced probability density, P n (x) — f£ p n (Qd£ for 
their distribution function and fJ, n {X) = J x p n for their measures. The iteration 
(2.4) induces a Frobenius- Perron (Chapman-Kolmogorov) equation for the distribution 
functions, 

J dh p{h)P n ^{A-\x - hj) = ]T ^P n _ x {f-\x)), (2.7) 

and accordingly for the densities and measures. The symbols f± denote the functions 
f±(x) := A(x) ± h. The Frobenius-Perron equation has a unique invariant measure 
fj,^ and the measures fin^ converge to in the weak topology of Borel measures 
on K. The invariant measure therefore is the measure of the effective field x in 
the thermodynamic limit b — > oo (n — > exo in the notation of the RIFS) . 

In this paper we focus our interest on the local magnetization in the bulk which 
is given by m^ k = (sj)jv at some site a < i < b inside the chain. To obtain (sj)jv we 
rewrite the partition function to a one-spin partition function with remaining spin Sj, 

Z N = £ exp(/3[(,f) + ^ ) ) Sj+ gs( y f))+ £ B(x^)}), (2.8) 

«i=±l j=a j=i+l 
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with two effective fields 

xf ] =A{xf + \) + h h i<j<b, x™=0, (2.9) 

=A{yf_\) + h h a<j<i, y { a N \ = (2.10) 

from the right and the left of site i respectively. The local magnetization at i is thus 
given by [4]- [6] 

= {si)N = tanh/3^ + A(yW)). (2.11) 

Let us introduce the notation 

f#M{X) := M-^X)) (2.12) 

for the mapping on Borcl measures induced by a measurable function /, e.g. 
tanh/3#(/i)(X) = /x(l//?artanh(X)). For the measure of m^ k we obtain in this 
notation 

M P=tanh/3 # (^ (a;) *A # ^)) (2.13) 

with I = i — a — 1 and r = b — i. As the effective fields share the same Frobenius-Perron 
equation (2.7) and the invariant measure of this equation is unique, the measures /i^ 
of the right-hand effective field in the thermodynamic limit b — > oo and /i^' of the 
left-hand effective field in the thermodynamic limit a — > — oo are identical. Therefore, 
as we will see below, the measure fif^ of the local magnetization in the bulk converges 
to 

// m > = tanh/? # ( M ( x > * A # ^) (2.14) 

in the thermodynamic limit a — > — oo, b — > oo (cf lemma 3.1 below), i.e. the 
local magnetization m^ k converges in distribution to a random variable m bulk with 

measure ^S m \ Please note that // m ) is space independent because of the uniqueness 
of the invariant measure of the Frobenius-Perron equation and the continuity of the 
convolution. 

The local magnetization at the boundary on the other hand is obtained if we 
consider only one effective field, i.e. 

m boundary = ^ = tanh ^ ( 2 .15) 

which has the measure tanh /3#fi^ . The multifractal properties of fj,^ are well known 
(cf [9, 11]) and general arguments show that tanh/3# has no effect on the Z? 9 -spectrum 
(cf [40]) such that the results apply to the measure of TO boundar y as we ll. The main 
point of this paper is the generalization to the magnetization of the bulk which is of 
greater physical interest. 



3. Convolution of multifractals 



In this section (/i„)„eN and (^n)neN denote sequences of bounded Borel measures on R 
which are Cauchy sequences with respect to the Hutchinson metric c?Hutch (cf [32] ) . As 
the space of bounded Borel measures on R is complete with respect to clutch [32], (fi n ) 
and (v n ) converge and we write fx := dHutch-limjwoo [i n and v := dHutch-linijwoc v n . 
As explained above we are interested in the properties of the convolution of bounded 
Borel measures. The convolution of two bounded Borel measures /i and v is always 
well defined (cf [41]) and will in the following be denoted by fi * v. As a first step we 
show that the convolution is continuous with respect to clutch- 



(3.3) 
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Lemma 3.1. Continuity of the convolution 

Let (mj)jgN an d (^»)»eN °e two monotonically growing unbounded sequences of natural 
numbers. Then (fj, m . * v n i)i£N converges to a bounded Borel measure in Hutchinson 
topology and the limit is clutch -linii—oo fJ, mi * v ni = u. * v. 

Proof. Let s > 0. The convergence of (fi n )neN and (v n )neN implies the existence of 
numbers M, N € N such that for all i > M clutch (Mm; , m) < e an d for all i > N 
duutchi^m, v) < e. Let N := max(M, N). For all i > N we then have 

sup< / f(z) 

(dz) - I f(z) pi * v(dz) | Lip(/) < 1 J . (3.1) 

The definition of the convolution of two measures implies J f(z)\i mi * v ni (dz) = 
JJ f(x + y)u. m .(dx)v ni (dy) and J f(z) pt* v(dz) = JJ f(x + y)u.(dx)v(dy). Inserting 
= - // f(x + y)/j,(dx)v ni (dy) + JJ f(x + y)v ni (dy)/j,(dx) we obtain 

duutchiUrm * V ni ,fl*l>) = 

sup < ( f(x + y)(J, m . (dx) - / f(x + y)n(dx) J v ni (dy) 

ff f \ \ (3 ' 2) 

+ J U fix + vWAdy) - J f(x + y)v{dy)\ u.(dx) Li P (/) < lj 

< J sup |y f(x + y)fi mi (dx) - j f(x + y)n(dx) Lip(/) < 1 j v ni (dy) 

+ y sup S^J f{x + y)u ni (dy) - J f(x + y)v(dy) | Lip(/) < lj p(dx). 
As the condition Lip(/) < 1 is translationally invariant we further obtain 
sup | J f(x + y)fi mi (dx) - J f(x + y)n(dx) Lip(/) < lj = 

sup | J f(x)fj, mi (dx) - J f(x)fj,(dx) Lip(/) < lj = d Hutch (n mi , n) < e. (3.4) 
In the same way 

sup^J f(x + y)v ni (dy) - J f (x + y)u(dy) Lip(/) < lj = d Hu tch(^n 4 , v) < £■ (3-5) 
We thus arrive at 

^Hutch^m; * v n ^n. * v) < J s u ni (dy) + J e fj,(dx) = (\\u ni \\ + WnW) s (3.6) 

in which \\i/ ni || = v ni (R) and = /z(R) denote the total mass of v ni and it respectively. 

□ 

As the metric clutch topology and the weak topology coincide on bounded Borel 
measures with compact support [32] we have the following corollary. 

Corollary 3.2. If supp fi and supp v are compact, w-linij-,^ \i m . * u n . = fi * v. 
Furthermore supp p*v is also compact. 

For the situation of the two-sided random-field Ising chain considered in this paper 
lemma 3.1 implies that the thermodynamic limit l,r —> oo can be carried out in an 
arbitrary way and that the result is the same as when first taking the thermodynamic 
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limit for the effective fields and then calculating the measure of the local magnetization. 
Having detailed knowledge of the properties of the D g -spectrum and the pointwise 
dimensions of the invariant measure of the effective field it is now interesting to 
investigate the relationship of the D 9 -spectra and pointwise dimensions of \i and v to 
the Dq-spectrum and pointwise dimensions of [i * v. The following lemmata allow us 
to transfer the knowledge about the multifractal properties of the invariant measure 
of the effective field gathered in [9] to the measure of the local magnetization. In 
the following we consider only measures with bounded, i.e. compact support. Let 
X- := min supp jj, > — oo denote the left and x+ := max supp ji < oo the right 
boundary of supp [i. For the boundaries of supp v we write y_ and y + . The pointwise 
dimension at the boundary of the support of /i * v can be obtained from the pointwise 
dimensions at x+, y + and y_. 

Lemma 3.3. Pointwise dimension of fi * v at the boundary of its support 

The left boundary of \i*v is z_ = X- + y- and the pointwise dimension D p (z-; = 
D p (x-]ii) + Dp(y^- 1 i i ). The result for the right boundary is analogous. 

Proof. The pointwise dimension of /i * v at z_ is defined as 

^(z_;^,) = lim ln(/i *^ (z - ))) (3.7) 
p e^o lne 

in which \i * v{B e (z-)) is given by 

V*v(B E (z-)) = J J l Bt {z-)(x + vMdx)v(dy). (3.8) 

The symbol lx denotes the characteristic function of a set X, i.e. lx{%) = 1 if x e X 
and = otherwise. The area in which l_B e (z_)( a; + y) is non-zero is shown in figure 1. 
Neglecting the regions in which either /i = or v = (or both), the relevant region is 
the dark gray triangle. As \i and v are positive measures, integration over the small 
square gives a lower and integration over the larger square an upper bound: 
rv-+i r x - + i ry-+ £ r x -+ e 

/ v{dy) ix(dx) <n*v{B e (z-j) < / v{dy) fi(dx). (3.9) 

J y- Jx- J V— Jx- 

Taking into account that /i = on (x- — s,X-) and v = on (y_ — £,y~) we can 
write 

v{B h ( y -))tJL{B h {x-)) <n*v{B e (z-)) <v{B e {yJ))n(B e (x-)) (3.10) 
to finally obtain 

In i/(5 t (?/_))+ In /iQB s (x_)) > hifi* u(B e (z-)) 
lnf+ln2 ~ lne 

> ]nv(B e (y-))+]ny(B e (x-)) 
lne 

which completes the proof as both sides of the inequality converge to D p (x_;/i) + 
D p {y_; v) as e — > 0. The proof for the right boundaries is obtained by applying the 
same arguments to fi(X) := /j,(—X) and v{X) := u(—X). □ 

To apply lemma 3.3 to the measure of the local magnetization in the ID RFIM 
it is important to know how the mappings A# and tanh/3# in (2.14) influence the 
pointwise dimensions of It turns out that they arc of no significance in this 

context because the pointwise dimension of the image measure at the image of some 
point is the pointwise dimension of the original measure at this point if the map under 
consideration is bi-Lipschitz. 



Convolution of multifractals and the local magnetization 



7 




Figure 1. Illustration of the proof of lemma 3.3. The diagonal strip is the region 
in which ^B e (z-){ x + v) ls non-zero. Therefore, the dark grey triangle is the 
relevant region with non-zero contributions to the integral (3.8). As fj, and u are 
positive measures, integration over the small square of side length | provides a 
lower and integration over the larger square of side length e an upper bound on 
the integral. 



Lemma 3.4. Stability of D p with respect to bi-Lipschitz maps 

Let f : M — > K be a bi-Lipschitz function and /i a bounded Borel measure on M. Then 
D P {f{x);f # {n)) = D p {x;u). 

Proof. As / is bi-Lipschitz so is / _1 and therefore 

L-'ly -x\< \f-\y) - f-\x)\ < L\y - x\ (3.12) 
for some constant L > 1. Then 

\f-\f(x) +e) - x\ = \r 1 (f(x)+e) - f-\f{x))\ < L\f(x)+e- f(x)\ = Le (3.13) 
and 

\x - r\f(x) - e)\ < L\f(x) - (f(x) -e)\= Le. (3.14) 

This implies f~ 1 {B e (f(x))) C Bl £ (x). In the same way one obtains B L -i e (x) C 
f-\B e {f{x))) such that 

lnn.(B L - lE (x)) > \nf #f i(B e (f(x))) > ln M (£? Le (*)) f3 15) 

The left and the right hand side of the inequality converge to D p (x; /i) such that the 
middle part which converges to Dp(f(x);f#/j,) also converges to this limit. □ 

For bounded measures with compact support it is sufficient that the function 
/ is bi-Lipschitz on an interval containing the support of the measure. As A(-) 
and tanh/3(-) are bi-Lipschitz on any finite interval lemma 3.3 and lemma 3.4 
directly imply L> p (m_; /i (m) ) = D p {m^; tanh /3 # (^ * -4#M (x) )) = D p {x-\^) + 
D p (A(x-):A # fi (x l) = 2D p (x-;n.^). In [9] lower (upper) bounds on D q for q < 
(q > 0) based on the pointwise dimension at arbitrary points in the support of the 
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Figure 2. Illustration of the main ideas of the proof of theorem 3.5. Figure 
a) applies to q > and figure b) to q < 0. The diagonal strip in a) represents 
the region of integration for (p, * v)i, the measure of one of the disjoint intervals 
of length e covering supp /i * v. The integral over the dark grey squares (and 
diagonally translated disjoint copies) provides a lower bound on (fi * v)i used in 
the case q > 1. Considering additionally the integral over the dashed squares 
gives an upper bound on (fi * v)i needed in the case < q < 1. 
The wide diagonal strip in b) is the region of integration for (n*v)i, the measure 
of one of the (intersecting) enlarged intervals of length 3e covering supp fi * v. 
The narrow dashed strip is the region of integration for the corresponding inner 
interval of size e. Integration over each of the six overlapping large squares of 
side length e (solid lines) and their disjoint by (— ne,ne) diagonally translated 
copies gives a lower bound on (p, * v)i such that the sum of the six integrals gives 
a lower bound on 6(fi * v)i. The narrow strip is contained in the union of all 
the interior small squares of side length | (dashed lines) assuring that the lower 
bound obtained is non-zero whenever the integral over the narrow strip is. This 
is an important point in the proof. The details are given in the text. 



measure were developed. These bounds can directly be applied to the Z? 9 -spectrum 
of /i^" 1 ' resulting in 

D q (^)>^L-2D p (x_)= q 21 * 2 (q<0). (3.16) 
q — 1 1 — qmA'(X-) 

This bound is a tight bound as long as the pointwise dimension at the boundary is 

weak. This is the case as long as D p {m ) > 1. The critical value h{ m '' 3 ^ determined 

by this condition is 

:1° L-i, „ 2flJ with (3.17) 



2/3 yR^ + e 2 ?^ 

R = 3 sinh(2/3J) - e~ 2/3J + \J (3 sinh(2/3 J) - e-^ iJ ) 2 - 1. (3.18) 

The critical value can also be interpreted in terms of the measure density. At this 
value of h the measure density at the boundary of the support changes from (for 
h < hi m ' 3) ) to oo (for h > h^ 3) ). 

For q > the corresponding bound is not of interest as the smoothness of in 
the region of small h implies smoothness of this region and thus D q = 1 for 
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all q > 0. The same applies to the connectedness of the support implying D = 1. 

The formerly discussed [9, 11] transition in the density of the effective field in 
which the slope of the coarse grained measure density at the boundary of the support 
changes from to oo also has an analogue. This effect occurs for the coarse grained 
measure density of the magnetization at £> p (m_) = | corresponding to 

h^** = -Lrsinh(2-I(l - 9e~ 40J )i) = . (3.19) 
P 

Note that the measure density of the local magnetization changes its slope at the 
boundary at the same critical value at which the measure density of the effective field 
at the boundary changes from to oo, cf [9, 11]. 

Apart from the relation between the pointwise dimensions of the convolution and 
its factors discussed so far there also exists a general relation between the D 9 -spectra. 

Theorem 3.5. Upper bound on D g (/x * v) 

The D q spectrum of the convolution is bounded from above by the sum of the D q - 
spectra of the factors, 

D q (n*v) < D q (n) + D q (v). (3.20) 

Proof. We need to distinguish three cases, q > 1, < q < 1 and q < 0. For the first 
two cases the improved multifractal formalism with enlarged boxes coincides with the 
usual one and for simplicity we will use the later in these cases. Throughout the proof 
sums of the form \i\ extend over alH g2 with > 0, i.e. boxes with zero measure 
are not taken into account (which is important for q < 0). Let e > 0. We denote 

Let q > 1. For any i e Z 



(fi * v)i := fi * v(Bz(x 2 i)) = J J ls § (x 2l )( x + y)fJ-(dx)u(dy) (3.21) 

is the integral over the diagonal strip in the (x,y)-plane shown in figure 2(a). 
Integration over the dark gray squares provides a lower bound on this integral. 

ifi*v)i > Y,^ B ^i+j)>i B ^i-o))- (3-22) 

3 

Taking the q-th power of both sides and using {J2i x i) q ^ ^2i x i for g > 1 and any 
positive numbers Xi we obtain 

(/**!/)? > 5>(£i(*2i +J -))MSi (X2i-j)y. (3.23) 

3 

For (fi * v)' q :— fi * v(B§. (x 2 i+i)) q we have in the same way 

(M*^)- 9 >J2^ B ^+3+i)) q ^ B ^2 l - 3 )) q - (3.24) 

3 

We denote fa := fi(B^(xi)) and i/j := v(B^(xj)). Summing (3.23) and (3.24) and 
over all i we get on the right hand side J^i J2j fcVj- ^ is straightforward to show that 
J2 i (fi * v)' q < 2 q+1 (cf Appendix A) such that the left hand side of the sum 

of (3.23) and (3.24) summed over all i is less than or equal to (2 q+1 + 1) ^2 Aft * v )1 ■ 
Taking the logarithm, dividing by lne and multiplying with l/(q — 1) we obtain 

1 lnE,(M*^+m2q+ 1 1 In ^ A + In g, v] 



In £ ~ q - 1 In f + In 2 
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which completes the proof for q > 1 as the left hand side converges to D q (fj, * u) and 
the right hand side to D q (fi) + D q {y) as e — > 0. 

Let < q < 1 and igZ. We again write (fj, * v),- L := fi * v(B^(x 2 i)), Hi := ^(-B|(xj)) 
and i/j := ^(£?| (a;.,-)). The solid and dashed squares in figure 2(a) and by (ne, —ne) 
diagonally translated disjoint copies cover the diagonal strip over which we need to 
integrate to obtain (/j, * v)i. We therefore have the upper bound 

l 

(u * v)% < ^2 ^2 V2i+j+kV2i-j- (3.26) 
j fe=-i 

Taking the q-th power and using (J2i x i) 9 — J2i x l for g < 1 and arbitrary positive 
numbers Xi yields 

l 

j k=-i 

When summing over all i each combination faVj appears at most twice in the sum on 
the right hand side such that 

< 2^ £>?i/?. (3.28) 

i i j 

Taking the logarithm of both sides, dividing by lne and multiplying with l/(q — 1) 
results in 

1 hi]T^(M * v )\ < 1 lnE^I+lnE J ^ + ln2 
q-l lne " g-1 In f +ln2 ' 1 ' ' 

The limit e -> yields D q (p * v) < D q (fj) + D q (v). 

Let q < 0. In this case we need the improved multifractal formalism with enlarged 
intervals. We use the notation 

By this choice we enlarge the e-intervals by e on both sides corresponding to k = 1 in 
Riedi's notation [40]. Furthermore we denote 

f MB.^)) (/i(B|(^)) > 0) f (a*)) MB.^)) > 0) 

' 4 ' \ (otherwise) ' 4 " - \ " (otherwise) ' { ' 

i.e. the |-intervals of [i and ^ are enlarged by | corresponding to k = \. This choice 
facilitates the proof and has no influence on the resulting D q as Riedi has shown (cf 
[40]). Let i G Z with (p * v)i > 0, i.e. the integral over the i-th interior interval 
is non-zero. When calculating (/j, * v)i we integrate over the wide diagonal strip in 
figure 2(b). The large squares B^(x 2 i+2j) x B^(x 2 i-2j), j € Z, are disjoint and are 
all contained in the strip. Therefore, 

(p * > ^2n2i+2jV2i-2j- (3.32) 

j 

This applies analogously to the other shown five squares and their by (ne, — ne) 
diagonally translated disjoint copies such that 

l 

6(p*v)i > X! M2i+j+fe^»-j- (3.33) 
j fe=-i 
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The integral over the narrow diagonal strip determines that (p * v)i is greater than 
zero. In the same way the integral over the small squares determines whether the 
terms on the right hand side are greater than zero. As the narrow strip is contained 
in the union of the small squares, the right hand side is greater than zero as (fi * v)i 
is. We therefore can take the g-th power on both sides and (omitting all terms being 
zero) use Xi) q < E, x\ for q < 1 and arbitrary positive numbers Xi to obtain 

6'(/* *«>)?<£ E Ai +i+ Ai-r (3-34) 
j fe=-i 

When summing over all i with (fi * v\ > 0, each combination /j^uj appears at most 
twice. Furthermore, adding terms which do not already appear only enlarges the right 
hand side. Therefore, 

6«5>*")i < 2 E^E4 (3-35) 

i i j 

From this we immediately obtain 

1 lnE^^+lne* 1 In + In + In 2 

q-1 Ins - q-1 ln| + ln2 1 ' ' 

which implies D q (fi * v) < D q (fi) + D q (v) in the limit e — > 0. □ 

Note that this proof easily generalizes to measures on W 1 . In [40] the invariance of 
the D 9 -spectrum with respect to bi-Lipschitz maps was shown, i.e. if / : M — > K is a 
bi-Lipschitz map then 

D q (n) = D,(/ # (/x)). (3.37) 

As in the case of the pointwise dimension it is sufficient that the function / is bi- 
Lipschitz on an interval containing the support of [i. Therefore, we can immediately 
deduce from theorem 3.5 and (3.37) that 

D q (^) = D q (^UA#^) 

< D q (/^) + D q (A # ti (x) ) - 2D q (p^). (3.38) 

As the -D g -spectrum of the invariant measure of the effective field is - at least on a 
numerical level - very well known (cf [9]) this provides interesting insights for the 
D g -spectrum of the measure of the local magnetization (cf figure 7) . In section 5 wc 
will discuss how to obtain the D g -spectrum of the measure of the local magnetization 
numerically and we will compare the results to the bounds obtained in this section. 



4. Convolution of measures on Cantor sets 

Letci 0) - [-\,\] andd n) be defined inductively by ci" } := MCi n ~ 1) )Uf a -{Ci n ~ 1) ) 
with f a +(x) = arr+irp and /„_ = ax— ^ Sl - The infinite intersection C a := n^Lo^"™^ 
is the a-Cantor set. On the approximating sets we define the probability densities 

Pi"» = fpfc^/a+W) + ^P^/a-V))- (4-1) 

The corresponding measures are denoted by fj^p(X) = J x p^ldx for any X G £>(K). 

The measures /zi™p converge to a limit measure \x a . p which is often referred to as an 
a-Cantor set with weights p and 1 — p. For a generic choice of a and p the measure 
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UL 




0.6 -o.: 




Figure 3. Illustration of the convolution of p^a} p with A#/iJ^p and p^a)p with 



p = j and A = =. The condition for the disjointness of 
the trapezoids, a/(l — 2a) < A < (1 — 2a) is clearly fulfilled for this choice of 
parameters such that in the limit n — > oo example 4.1 applies. 



fj, a ,p is a multifractal. For an illustration cf figure 3 and 4. In the following example 
we calculate the D q spectrum of the convolution of fj, a ,p with A#^ a ,p, a "compressed" 
version of itself (A < 1). As this is in general a hard problem we discuss two examples. 



Example 4-1. Let < p < 1 and a/(l — 2a) < A < (1 — 2a). This is meaningful 

for a < j. Then, the two intervals of supp ^a)p fit into the gap of supp \#/i^ p 

and on the other hand the complete supp \#)ia,p fits into the gap of supp (J-^p. 
Self-similarity of C a and \#C a imply that for any given n and y £ R at most 

one pair of bars in pi™p(x) and A#pi"p(y — x) can overlap. The convolution is 
therefore a collection of trapezoids as shown in figure 3. We denote the intervals 
of supp fia^p * A#yui"p (the bases of the trapezoids) by T^ n \ j — 1, ...,2". For 
larger n only the structure of the measure within the trapezoids but not their total 
measure changes, i.e. ^i a ,p * A#/Lt 0iP (T- n ^) = * X^^l(T^). We therefore can 

use [i a ,p * A^/x a .p and fj^fy * A#/ii"p interchangeably. This fortunate circumstance is 
due to the fact that the self-similarity of the Cantor sets induces a direct iteration 
for the convolution making it self-similar itself. The analytical treatment of the D q - 
spectrum in this example is essentially based on this fact. If we choose e„ := a" + Aa", 
which is the width of the trapezoids at level n, boxes B3 Sn (xi), Xi — is n , i S Z with 
fi a , P (B£s, (xi)) > contain at least one whole trapezoid of level n and intersect at most 
four. Thus, denoting 

Va,p * A #AVp(-B§e„ {Xi)) (fJ, atP * \ # ll a ,p(B^ (x t )) > 0) 

(otherwise) 
we obtain for q > 0, q ^ 1 



(4.2) 



fjb a ,p * A#Mo,p(7 1 . 



< H\ < (4 max Ma , p * X^ a . p (T^))\ (4.3) 

where is the index of a trapezoid completely contained in -B| £n (x^ and J(i) is the 
set of the indices of all trapezoids intersecting Bz € (x^. As any trapezoid can appear 
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Figure 4. Illustration of the convolution of fi~a)p with itself and fi~a} p with itself 
for a = j and p = |. For this choice of parameters example 4.2 applies in the 
limit n — > oo. 



at most four times on the right hand side when summing over i and any trapezoid 
appears at least once on the left hand side this implies 

E * ^AT^Y <J2^< 4E( 4 ^ * ^a, P {^ n) )) q (4.4) 

3 * 3 

The measures of the trapezoids can explicitly be calculated such that 

E * *#*>AT} n) ) q = E (?) f?) (P^ 1 - P)"" V(l - P) n - l ) q - (4-5) 



fc,J=0 



Applying £ fc (fe)(p 9 ) fe ((l -p)T~ fc = (p 9 + (1 -p)T we obtain 

(p 9 + (l-p)«) 2n <E^i <4«+ 1 (p« + (l-rt") 2 " (4.6) 

and therefore 

For q < the argument is the same with reversed inequality signs which leads to the 
same result. 

For q = 1 we calculate the limit q — ► 1 of (4.7) yielding 

5— >i q — 1 In a 

= 2(pLnp+(l -p)ln(l-p))/lna. (4.8) 

For any A G [a k+1 /(l — 2a) 1 a k (l — 2a)], k G N, the arguments above apply to 
Ha^p^ * A#/xi™p which according to lemma 3.1 also converges to \i a , v * A#/i aiP . 
Therefore, (4.7) and (4.8) apply to all A taken from these intervals. For an example 
cf figure 5. 
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Example 4-2. Let 0<p<l,a<| and A = 1. We then choose e n := 2a™ and boxes 

of length s n in such a way that each box covers one of the spikes of /zi^p * A*o?p- This 
is a permitted choice and avoids any complications with q < such that we refrain 
from using enlarged boxes here. The situation is like shown in figure 4, i.e. for n = 1, 
* Ha~}p)i = (1 — p) 2 , 2p(l—p) and p 2 for i = 1,2 and 3 respectively. As in example 
4.1 self-similarity implies for our specific choice of boxes that we can use fi^p * ji^p 
and Ha iP * ^ aiP interchangeably as they coincide on all boxes. Furthermore, the result 
of the convolution of the next iteration, n = 2, can be constructed by replacing each 
triangle in figure 4 by the complete figure and choosing the corresponding weight. 
Therefore, the (/J,a}p * Ha}p)i, i = 1, . . . , 9, sum up to 

E E (f) ( 2 1 k ) o- p) 2k (Mi p)r k p 2(2 - k - l) (4.9) 

fc=o ;=o ^ ' ^ ' 

= {(l-p) 2 + 2p(l-p)+p 2 ) 2 . (4.10) 
This continues for larger n such that we have 

n n—k / \ / j\ 

E«*^ = EEu)( n 7 ) ((i- P ) 2fc (2 P (i-rtr-v ( - fe - i) ) 9 (4.11) 

= ((1 - p) 29 + 2V(1 - P) 9 + P 29 )"- (4.12) 
For g 7^ 1 this yields 

D = _J_ lim »ln((l-p) 2g + 2V(l-p)^+ P 2 g) ^ 
9 g — 1 rt->oo n In a + In 2 

1 k((l-p) 2 » + 2^(l-p)Hp 2 ') f4U1 
o-l lna ' 1 ' J 

The limit q — > 1 results in 

.. 1 ln((l-p) 2 9 + 2V(l-p) 9 +P 29 ) 
L)\ = lim 

9-»i q — 1 lna 
= 2(plnp + p(l -p)ln2 + (1 -p)ln(l -p))/lna. (4.15) 

As in example 4.1 the arguments can be repeated for the a"-fold value of A such that 
(4.14) and (4.15) are also correct for all A = a™, n <E N. This means that we can 
calculate the D 9 -spectrum at specific points between the intervals where it is given by 
(4.7) and (4.8). The results are illustrated for a — | and p = | in figure 5. 

In the case of the random-field Ising model we are interested in the dependence of 
the -Dg-spectrum on the strength h of the random field which rather corresponds to 
varying a in the convolution of Cantor sets. Viewing (4.7) and (4.8) in this light we 
can choose A = ^ such that for all a < j example 4.1 applies. This results in the 
right part of the Z?,j-spectrum shown in figure 6. For the left part with a close to one 
we have the usual lower bounds for q < based on the pointwise dimension at the 
boundary of the support and D q = 1 for q > 1 from the regularity of fJ, a ,p- 

The convolution of two-scale Cantor sets shows a far richer behaviour than the 
two examples above. It turns out that the lacking strict self-similarity does not allow 
the kind of analytical treatment we used so far. Numerical investigations show that 
the D 9 -spectrum strongly depends on the two scales of the Cantor sets. Furthermore, 
in a large parameter regime the obtained numerical estimates are extremely sensitive 
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Figure 5. D q spectrum of fi a ^ p * \#n a ,p as a function of In A for q = —20, —6, 
—3, —2, —1, 0, 1, 2, 4 and 20. The thick lines are the exact D q obtained in 
example 4.1 and the points the D q obtained in example 4.2. The thin lines are 
numerical results obtained from iteration depths I = r = 8 compared to I = r = 7 
in the new natural partition method (cf section 5.3). 




Figure 6. Dq-spectrum of /i a ,p * A#/i aiP as a function of —In a with P = §, 
A = \ and q = -20, -6, -3, -2, -1, 0, 1, 2, 4, 20. The lines to the left of 
— In a = In 4 are the usual lower bounds based on the pointwise dimension at the 
left boundary of C a ,p- The lines on the right arc the calculated exact values of D q 
according to (4.7) and (4.8). The points are numerical results for iteration depths 
I = r = 6 compared to I = r = 5. 
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to the iteration depth to which the Cantor sets are generated. The situation is of 
similar complexity as for the measure of the local magnetization in the one-dimensional 
random-field Ising model discussed below. We therefore refrain from discussing this 
case and turn to the characterization of the measure of the local magnetization. 

5. D q -spectrum of the local magnetization 

In this section we present our numerical methods and the resulting D g -spectra for the 
measure of the local magnetization in the random-field Ising model. We obtain the 
D g -spectrum of the measure of the local magnetization in several steps. 

5.1. Generation of the measure of the effective field 

The measure of the local magnetization in the bulk of the random-field Ising chain 
is given by (2.11), i.e. it is the convolution of the measure of the effective field 
x n with a distorted version of itself. Therefore, as a first step for a numerical 
treatment, an approximation of this measure is needed. It turns out that a very 
useful approximation is the following (cf [9, 11]). We take the partition of the 
invariant interval I = [a;* given by the points {xi} = {f{ a }„(x*_), /{ CT }„(x+)}. 
It is sometimes called the "new natural partition" (cf [11]). The points x*_ and x* + 
are the fixed points of the iteration (2.4). The symbols f{ a } n ( x ) denote the iterated 
functions f {a}n (x) := f ai o f a2 o . . . o / CTn (x) with f + (x) = A(x) + h, /_ (x) = A(x) - h 
and {a} = {<7i, 02, . . . , cr„} with <7j G {+, — }. The n-fold application of the Frobenius- 
Perron equation (2.7) on some initial distribution P Q yields [11] 



We then calculate the measure /x„ ([xj, x^+i]) = P n {x i+1 ) — P n {xi) on the intervals of 

(x) 

the new natural partition and approximate the density p„ ' as constant within these 
intervals. The resulting histogram is used for the convolution. 

5. 2. Convolution 

The next step toward the D g -spectrum of the measure of the local magnetization is to 
calculate the convolution (2.14). We will need to calculate the measure of some given 
intervals, say [m„m i+ i]. It is given by 




(5.1) 




^iy([rni,m i+ i]) = tanh /3 # (nf * A # ^ x) )([m i; m i+1 ]) 



l//3tanh _1 (m i+ i) 




(5.3) 



1//3 tanh" 1 (mi) 
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For and p\ x ' we substitute our piecewise constant approximations pf and pf ' and 
denote the inner integral as F(y). This function is piecewise of the form KiA(y) + rji 
and thus can be represented by th coefficients Ki, rn and the endpoints yt of the 
intervals [yi,yi+i] on which the particular Ki and r]i are valid. The outer integral in 
(5.3) is therefore approximately 

£ r +1 ^ x \y)(iuA{y)+rH)dy. (5.4) 

As is piecewise constant this integrals can easily be calculated provided J A(y)dy 
is known. This integral is given by 

J X A(y)dy = JL^-e 2 ^- 7 )) - Li 2 (-e 2 ^+ J ))) -Jx + C (5.5) 

in which Li 2 denotes the second polylogarithmic function:):. It turns out however 
that for the short intervals we need to integrate on the implementation of the 
polylogarithmic function Li 2 for double precision numbers is less precise than a simple 
fifth order Taylor expansion of / A(y)dy around the center of the intervals [yi,yi+i]. 
We therefore use the expansion in our numerical studies. 

By the methods described thus far we have gained the ability to obtain the 
measure [iff of any given interval [mj,TOj+i] and not too large iteration depths I 

and r. We now use two different methods to estimate the -D g -spectrum of [i^ based 
on this. 



5.3. Determination of the D q - spectrum 

The first method uses coverings of supp [iff with boxes of equal size. We choose 
boxes of size et ■— £o • s k , k — 1, 2, . . . , N, for some s < 1 and points Xi — isk, i € Z. 
We then calculate 

Z[ B) =Y.^\B hk {x i ))Y (5.6) 

W 

in which only indices i fulfilling [iff (B^(xi) > are considered. A linear fit of 

\nZ k ' as a function of lne^ yields hiZ^, ; ~ r q \nek + const, providing the desired 
estimate D q = r q /(q — 1) of the generalized fractal dimension D q . 

The second method is based on the stationarity of a suitably chosen partition 

(x) (x) 

function. When observing the process of the convolution of [i\ and A#[T r more 

(x) 

closely it becomes clear that there is a qualitative change whenever bands of [i\ 

(x) 

and bands of A#[ir start or cease to overlap, i.e. for values m, of the magnetization 
obeying the condition 

-artanh(TOi) - g(xij) = x r ^ k (5.7) 

(x) (x) 

where xij and x r ^ are points of the new natural partition of [i\ and [if- respectively. 
This condition leads to 

m t =tanh(3(f {a}l (x* ± )+A(f {a]r (x* ± )))- (5-8) 

t The 2-nd polylogarithmic function is Li 2 (z) := J2T=l W = ~ Jo dt. 
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We will employ these points as a new natural partition for the measure of the local 
magnetization. It turns out that there exists a natural degeneracy within the set {m^} 
induced by the trivial identity 

A(a) + f-(b) = A(b) + f-(a) (5.9) 

and other such identities comprising higher iterations of /+ and /_ . These degeneracies 
have to be removed "by hand" by the algorithm. In the spirit of [42] we then define 
the partition function 

— / \ (5.10) 

(m i+ i - mi) T i 

on this new natural partition and determine D q from the condition 

\nZ hr (q, Tq ) -\uZ Vtr ,{q,T q ) = (5.11) 

with some iteration depths I, r, I 1 , r'. 

Figure 7 shows a summary of the obtained numerical results. In region a) the 
behaviour of the D q with negative q is dominated by the pointwise dimension at the 
boundary of the support of (j,( m ^ such that the lower bounds based on it coincide with 
the obtained numerical values. Furthermore, the numerical results are very stable for 
all iteration depths. We do not show the result of the box method in this region as it 
is known that box methods systematically underestimate D q for q < if it strongly 
depends on only few points in the support. The generalized dimensions for q > are 
all 1 because the measure is smooth. The numerical results are in perfect agreement 
with this statement. 

In region b) all D q for q < arc 1 for h > fti™' 3 ' '■ At some point they again arc 
greater than 1. In this region the method based on the new natural partition yields 
completely different results for different iteration depths. We therefore are not able to 
deduce the asymptotic behaviour from the scaling in finite iteration depth. Most data 
points had to be left out for this reason. Provided h is large enough (h > 1.7) the 
numerical results of the new natural partition method are again stable for all iteration 
depths which leads to small error bars in figure 7. For q > we have in regions b) and 
c) perfect agreement with the upper bounds obtained from the pointwise dimension 
at the boundary of the support of ^S m \ 

The difficulties in obtaining the asymptotic scaling for q < in the region 
1 ~ h ^ 1-7 are of the same type as encountered in the convolution of two-scale 
Cantor sets. This shows that this is an effect of more than one relevant scale present 
(infinitely many in this case). As the asymptotic scaling seems not to be attainable we 
have the impression that from a physicists point of view we should pose the question 
what an experimentalist would observe. In any experiment the scale of resolution is 
bounded from below. This corresponds to the situation of the box method where the 
scale is bounded by the size of the smallest box (whereas the scale in the new natural 
partition can get more or less arbitrarily small already for finite iteration depths). We 
therefore surmise that the D q estimates based on the box method are the physical 
results in this region. As it turns out the results of the box method are fairly robust 
against changes in the iteration depth whereas they depend on an appropriate choice 
of box sizes. As a rule of thumb the smallest box size should be of the order of the 
length of the longest band of the new natural partition. The results of the box method 
are shown as lines in regions b) and c) in figure 7. 



Zi, r (q,T q ) = V 




Figure 7. Numerical results for the D q -spcctrum of the measure of the local 
magnetization in the bulk. We considered q = —20, —6, —3, —2, —1, 0, 1, 2, 4 
and 20. All points with errorbars were obtained by the method based on the new 
natural partition. In regions a) and b) and q ^ the results of all combinations 
of iteration depths I = r = 8, I = r = 9 and I = r = 10 were used and the 
number representation was C++ long doubles. In region c) and also q ^ all 
combinations of iteration depths l = r = 5, l = r = 6 and I = r = 7 were 
used with numbers of the arbitrary precision library "CLN" with guaranteed 50 
decimal digits. The errorbars are obtained from the standard deviations of the 
average of the results of the three possible combinations of iteration depths. All 
points with standard deviation greater than 0.05 are not shown. For q = we 
used iteration depths up to I = r = 13. The dashed lines in a) are lower bounds 
on D q and the dashed lines in b) and c) are upper bounds both based on the 
pointwise dimension at the boundary of the support of /i' m ) . The dotted lines 
arc the twofold numerical results for the D q (q < 0) of the effective field (cf [9]) 
and are thus upper bounds on the D q (q < 0). Finally, the other lines in b) and 
c) are the results for D q (q < 0) obtained by the box scaling approach. The usual 
spacing in h for all data points is 0.02 except for the region between h = 0.56 
and h = 0.66 where we chose a finer spacing of Ah = 0.005 in the box method. 
P = J = 1. 



The local minimum of the D q , q < at h ~ 1.6 can be understood as a change 
in the overlap structure of bands of the new natural partition of the two convoluted 
measures for different h. For h « 1.5 there is a relative position for the approximations 
/Li[ x ^ and A^fii^ for which only the two very weak bands around and A(x* + _) 
overlap leading to a very weak band in the convolution, cf figure 8 for an example. 
This results in the large values of D q we observe numerically for the iteration depths 
I = r < 10. For h w 1.6 the band structure is of a form that no such position can be 
found in the iteration depths under consideration. For any relative position of the two 
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Figure 8. Illustration of the situation leading to the weakest band in the 

(x) (x) 

convolution of /i ; and A^fi r at h = 1.54 and I = r = 6. The integration 
over the product of (the density of /j,i x \x), upper part) and A^p^ (y — x) 

(x) 

(the density of A#p). , lower part) yields the density of the convolution at y. 
In this figure y = 0.309965 which is the position of the weakest band. Only 

for the weak band around x^ ^ the two densities are simultaneously non-zero 

(see inlay). This leads to a very small value for the density of the convolution. 
Note that even though the two convoluted measures are already rather sparse, 
the convolution still has non-fractal support at this h. (3 = J = 1. 



approximations of the measures more than one pair of bands or considerably stronger 
ones overlap. This results in the considerably smaller values of D q . For larger h the 
formation of a very weak band in the convolution reappears and we again get large 
values of D q {q < 0). For other iteration depths the situation can again change as 
there is no strict self-similarity of the measure fi^ m ' . From another point of view for a 
given iteration depth the values of D q strongly depend on the random field strength 
h due to similar changes in the overlap structure as discussed so far. We expect that 
on any iteration depth this situation is qualitatively the same. 

6. Conclusions 

In order to investigate the multifractal properties of the measure of the local 
magnetization in the random-field Ising chain we studied the relation between 
the multifractal properties of two measures and their convolution. The pointwise 
dimension at the boundary of the support of the convolution turned out to be the 
sum of the pointwise dimensions at the boundary of the two measures. This enabled 
us to calculate the pointwise dimension at the boundary of the support of the measure 
of the local magnetization and employ this to give bounds on the D g -spectrum. 

We furthermore were able to prove that the generalized dimensions D q of the 
convolution are bounded from above by the sum of the D q of the two convoluted 
measures. This yields upper bounds on the D g -spectrum of the measure of the local 
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Figure 9. Measure densities of the effective field (left column) for iteration depth 
16 and the local magnetization (right column) for iteration depth I = r = 8. The 
details are explained in the text. (3 = J = 1. 
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magnetization. 

The general results were illustrated for the convolution of Cantor sets with weights 
and our main application, the measure of the local magnetization in the random-field 
Ising model. We also performed numerical studies of the D^-spectra employing a 
box method and a method based on a new natural partition. The numerical data 
is consistent with the exact inequalities but reveals a complicated structure for the 
D 9 -spectrum of the measure of the local magnetization for q < and 1.2 < h < 1.7. 
The numerical instability of the method based on the new natural partition in this 
region can be understood by analysing the band structure of the approximations of the 
invariant measure of the effective field. Nevertheless, it prevents us from obtaining the 
asymptotic scaling. We therefore took a pragmatic approach and investigated what 
an experimentalist could possibly hope to observe resulting in the shown D q estimates 
based on the box method. 

Investigating the probability measure of the local magnetization of the one- 
dimensional random-field Ising model we for the first time considered a physical 
quantity which in principle is measurable. Furthermore, the general results about 
convolutions of multifractal measures can be of interest in other areas as well as it is 
not unusual that sums of random variables appear. 

All figures were presented for a generic choice of parameters, (3 = J = 1. The 
behaviour for f3, J — > and /3, J — > oo is more or less trivial. For (3 — ► the probability 
measure of the local magnetization is a Dirac measure at 0. For J — > it is the sum 
of two Dirac measures at ±tanh/3ft. In the limit (i — > oo the distribution of the 
effective field is a finite sum of Dirac measures [5, 6] and therefore the measure of the 
local magnetization is also such a sum. Finally, in the unphysical limit J — > oo the 
function A{x) is the identity, A(x) = x, and the RIFS for the effective field therefore 
ceases to be contractive. The measures \i n of the effective field for finite system sizes 
need not converge in Hutchinson metric in the thermodynamic limit n — > oo. The 
proof for existence and uniqueness of the convolution of lemma 3.1 does therefore 
not apply. One needs to consider finite systems and take the thermodynamic limit 
in the end. The RIFS is a symmetric random walk with step size h and the local 
magnetization is m^ r — tanh/3(xj + y r ) where xi and y r are two random walks of 
length I and r respectively. Thus, l//3artanhm; !? , is also a random walk of length 
n = I + r. One therefore can deduce that the probability for the local magnetization 
to take values in any closed interval X C [—1,1] tends to zero in the thermodynamic 
limit n — > oo because this corresponds to the probability for the random walk to stay 
in a finite region. The measure of the local magnetization thus converges to the sum 
of two Dirac measures at ±1 in the weak topology of Borel measures. In all cases no 
multifractal effects of interest can be observed. 

As a byproduct of our numerical algorithm we can easily produce the measure 
density of the local magnetization itself, cf figure 9. In the figure we show the 
measure density p^ of the effective field and the measure density p( m ) of the local 
magnetization for some values of the random field strength h. One can clearly see 
that p( m ) is much smoother than p^ i n accordance with the general belief. For 
h = 0.02 both measures are smooth and the slope at the boundary is zero. For 
h = 0.2 > the slope of p^ is already infinite whereas the slope of p^ remains 
zero. For h — 0.4 > = hi" 1 '^ the density p^ is infinite at the boundary whereas 
p( m ^ is zero but has infinite slope. For h = 0.7 > /i£ m ' 3 ^ the density p( m ^ also is infinite 
at the boundary. The fractality of the support is in the same way "delayed" for p( m ) . 
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For h = 1.0 > hi 1 ' the support of p( x ) is already fractal but the support of p 1 -" 1 ' is still 
Euclidean. 

Overall, there is a gradual transition from a monomodal strongly peaked 
distribution for small random field to an even more strongly peaked bimodal 
distribution for large random field. The local magnetization (which is a thermo- 
dynamic average but still a random variable with respect to the probability space of 
the random field) shows a transition from a paramagnetic situation where the most 
probable value is zero to a ferromagnetic situation where the most probable value is 
±1. Between these extremal situations lies the multifractal regime. The distribution 
always remains symmetric such that this is not a phase transition; there is no symmetry 
breaking even if a small homogeneous field is applied. 



Appendix A. 

Let Xi = ie, I G Z and x\ = Xi + y, i £ Z be two grids which are shifted by y with 
respect to each other and let q > 0. In this appendix we show that 

J2KB § W) q <2 q+1 J2v(B i &)) q - (A.I) 

i i 

Let i G Z and denote /Uj := /j,(B^(xij) and /U- := p(_B|(x-)). Clearly, < J2je.7(i) 
where J(i) — {j £ Z : B^(xj) n B^(x' i ) ^ 0}. As B^(x' i ) intersects at most two 
B*(xj) the set J(i) has at most two elements and we can write 

tfKVmaxtf. (A.2) 
jeJ(i) J 

On the other hand each Bz(xj) intersect at most two B^(xi) for a fixed j such that 
appears at most twice when summing (A.2) over all i. Therefore, 

5>r<2.2«5>« (A.3) 

i i 

which is the result claimed in (A.l). 
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